Search results

An elasto‐viscoplastic constitutive model valid under general complex cyclic loadings is proposed to describe the multiaxial behaviour of engineering metallic materials. Its numerical implementation in the framework of the finite element method is underlined. Its application to the prediction of some various tests is compared to the experimental responses.

Complex constitutive equations including non‐linear kinematic as well as isotropic hardening rules are introduced in the framework of the finite element method to analyse the behaviour of structures in quasi‐static elasto‐viscoplasticity under cyclic loadings. An implicit time marching process is developed together with a substructuring technique and the corresponding numerical details are underlined. Practical examples are treated to demonstrate the capabilities of the proposed finite element code.

The present paper aims to explore a computational homogenisation procedure to investigate the full geometric representation of yield surfaces for isotropic porous ductile media. The effects of cell morphology and imposed boundary conditions are assessed. The sensitivity of the yield surfaces to the Lode angle is also investigated in detail.

The microscale of the material is modelled by the concept of Representative Volume Element (RVE) or unit cell, which is numerically simulated through three-dimensional finite element analyses. Numerous loading conditions are considered to create complete yield surfaces encompassing high, intermediate and low triaxialities. The influence of cell morphology on the yield surfaces is assessed considering a spherical cell with spherical void and a cubic RVE with spherical void, both under uniform strain boundary condition. The use of spherical cell is interesting as preferential directions in the effective behaviour are avoided. The periodic boundary condition, which favours strain localization, is imposed on the cubic RVE to compare the results. Small strains are assumed and the cell matrix is considered as a perfect elasto-plastic material following the von Mises yield criterion.

Different morphologies for the cell imply in different yield conditions for the same load situations. The yield surfaces in correspondence to periodic boundary condition show significant differences compared to those obtained by imposing uniform strain boundary condition. The stress Lode angle has a strong influence on the geometry of the yield surfaces considering low and intermediate triaxialities.

The exhaustive computational study of the effects of cell morphologies and imposed boundary conditions fills a gap in the full representation of the flow surfaces. The homogenisation-based strategy allows us to further investigate the influence of the Lode angle on the yield surfaces.

This work addresses the computational aspects of a model for elastoplastic damage at finite strains. The model is a modification of a previously established model for large strain elastoplasticity described by Perić et al. which is here extended to include isotropic damage and kinematic hardening. Within the computational scheme, the constitutive equations are numerically integrated by an algorithm based on operator split methodology (elastic predictor—plastic corrector). The Newton—Raphson method is used to solve the discretized evolution equations in the plastic corrector stage. A numerical assessment of accuracy and stability of the integration algorithm is carried out based on iso‐error maps. To improve the stability of the local N—R scheme, the standard elastic predictor is replaced by improvedinitial estimates ensuring convergence for large increments. Several possibilities are explored and their effect on the stability of the N—R scheme is investigated. The finite element method is used in the approximation of the incremental equilibrium problem and the resulting equations are solved by the standard Newton—Raphson procedure. Two numerical examples are presented. The results are compared with those obtained by the original elastoplastic model.

The purpose of this paper is to systematically investigate the dynamic strain aging (DSA) effect in solid solution treated IN718 at different temperatures through experiments and simulations to gain an understanding of the inelastic deformation mechanisms.

In the present work, uniaxial tensile tests have been carried out in conjunction with finite element (FE) simulations to investigate the behaviour of the solid solution treated Inconel 718 superalloy at different temperatures and strain rates. Dynamic strain aging (DSA) effects, which manifested during the tests in the form of a negative strain rate sensitivity and stress serrations, are investigated. The most significant DSA effect occurs at 500°C and at a strain rate of 10–4 s-1. In a newly proposed rate-dependent constitutive formulation, the DSA model, proposed by McCormick, Kubin and Estrin, was introduced into slip-assisted solute hardening, and an activation energy-dependent exponential flow rule was adopted.

The observed negative strain rate sensitivity and stress serrations are well predicted by a 3 D FE. The FE results indicate that the equivalent plastic strain rate distribution in the specimen gauge length is as highly inhomogeneous as in the other materials exhibiting DSA effects such as aluminium and titanium alloy. During inelastic deformation, propagating high strain rate bands can be closely correlated to the stress serrations.

For the DSA effect in solid solution treated IN718, the existing researching mainly focuses on the mechanical properties experiment and microstructure observation. In this study, a constitutive formulation, combined with the DSA model, has been proposed, and the mechanical behaviors, including the DSA effect, have been well predicted by a finite element model.

The purpose of this paper is to suggest a robust hybrid method for updating the stress and plastic internal variables in plasticity considering damage mechanics.

By benefiting the properties of the well-known explicit and implicit integrations, a new mixed method is derived. In fact, the advantages of the mentioned techniques are used to achieve an efficient integration.

The numerical studies demonstrate the high precision and robustness of the suggested algorithm.

The perfect von-Mises plasticity together with Lemaitre damage model is considered within the realm of small deformations.

Updating stress and plastic internal variables are of utmost importance in elastoplastic analyses of structures. The accuracy and efficiency of stress-updating methods significantly affect the final outcomes of nonlinear analyses.

The idea which is used to derive the hybrid method leads to an efficient integration method for updating the constitutive equations of the damage mechanics.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from excessive mesh dependence when strain‐softening models are used in numerical analyses and cannot reproduce the size effect commonly observed in quasi‐brittle failure. In this contribution three different approaches will be scrutinized which may be used to remedy these two intimately related deficiencies of the classical theory, namely (i) the addition of higher‐order deformation gradients, (ii) the use of micropolar continuum models, and (iii) the addition of rate dependence. By means of a number of numerical simulations it will be investigated under which conditions these enriched continuum theories permit localization of deformation without losing ellipticity for static problems and hyperbolicity for dynamic problems. For the latter class of problems the crucial role of dispersion in wave propagation in strain‐softening media will also be highlighted.

Simple but also accurate models are needed to predict the failure response of concrete structures. Simplicity involves modelling assumptions while accuracy involves objectivity of both the experimentally identified model parameters and the numerica results. For concrete‐like heterogeneous and brittle materials, the modelling assumptions idealizing the material as a homogeneous continuum with classical linear or non‐linear behaviour, leads to some problems at the identification stage, namely the size effect phenomena. A continuum damage model, representing the non‐linear behaviour due to microcracking, is proposed here for predictive computations of structural responses. A Weibull based theory is used to determine, in a statistical sense, the value of the initial damage threshold. The essential influence of material heterogeneity on the damage evolution, is accounted for with a bi‐scale approach which is based on the idea of the non‐local continuum with local strain. It has already established that the non‐local approaches yield realistic failure predictions and the numerical results are convergent for subsequent mesh refinements. The applications presented here show the ability of the approach to predict the failure response of concrete structures without being obscured by size effect problems.

It has been usual to prefer an enrichment pattern independent of the mesh when applying singular functions in the Generalized/eXtended finite element method (G/XFEM). This choice, when modeling crack tip singularities through extrinsic enrichment, has been understood as the only way to surpass the typical poor convergence rate obtained with the finite element method (FEM), on uniform or quasi-uniform meshes conforming to the crack. Then, the purpose of this study is to revisit the topological enrichment strategy in the light of a higher-order continuity obtained with a smooth partition of unity (PoU). Aiming to verify the smoothness' impacts on the blending phenomenon, a series of numerical experiments is conceived to compare the two GFEM versions: the conventional one, based on piecewise continuous PoU's, and another which considers PoU's with high-regularity.

The stress approximations right at the crack tip vicinity are qualified by focusing on crack severity parameters. For this purpose, the material forces method originated from the configurational mechanics is employed. Some attempts to improve solution using different polynomial enrichment schemes, besides the singular one, are discussed aiming to verify the transition/blending effects. A classical two-dimensional problem of the linear elastic fracture mechanics (LEFM) is solved, considering the pure mode I and the mixed-mode loadings.

The results reveal that, in the presence of smooth PoU's, the topological enrichment can still be considered as a suitable strategy for extrinsic enrichment. First, because such an enrichment pattern still can treat the crack independently of the mesh and deliver some advantage in terms of convergence rates, under certain conditions, when compared to the conventional FEM. Second, because the topological pattern demands fewer degrees of freedom and impacts conditioning less than the geometrical strategy.

Several outputs are presented, considering estimations for the J–integral and the angle of probable crack advance, this last computed from two different strategies to monitoring blending/transition effects, besides some comments about conditioning. Both h- and p-behaviors are displayed to allow a discussion from different points of view concerning the topological enrichment in smooth GFEM.

The purpose of this paper is to discuss the recently developed multi-physics computational model for the conventional Gas Metal Arc Welding (GMAW) joining process that has been upgraded with respect to its predictive capabilities regarding the spatial distribution of the mechanical properties controlling the ballistic limit (i.e. penetration resistance) of the weld.

The original model consists of five modules, each dedicated to handling a specific aspect of the GMAW process, i.e.: electro-dynamics of the welding-gun; radiation-/convection-controlled heat transfer from the electric arc to the workpiece and mass transfer from the filler-metal consumable electrode to the weld; prediction of the temporal evolution and the spatial distribution of thermal and mechanical fields within the weld region during the GMAW joining process; the resulting temporal evolution and spatial distribution of the material microstructure throughout the weld region; and spatial distribution of the as-welded material mechanical properties. The model is upgraded through the introduction of the sixth module in the present work in recognition of the fact that in thick steel GMAW weldments, the overall ballistic performance of the armor may become controlled by the (often inferior) ballistic limits of its weld (fusion and heat-affected) zones.

The upgraded GMAW process model is next applied to the case of butt-welding of MIL A46100 (a prototypical high-hardness armor-grade martensitic steel) workpieces using filler-metal electrodes made of the same material. The predictions of the upgraded GMAW process model pertaining to the spatial distribution of the material microstructure and ballistic-limit-controlling mechanical properties within the MIL A46100 butt-weld are found to be consistent with general expectations and prior observations.

To the authors’ knowledge, the present work is the first reported attempt to establish, using computational modeling, functional relationships between the GMAW process parameters and the mechanical properties controlling the ballistic limit of the resulting weld.